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Abstract 

A new approach to quantum field theory at finite temperature and density in arbitrary space-time 
dimension D is developed. We focus mainly on relativistic theories, but the approach applies to 
non-relativistic ones as well. In this quasi-particle re-summation, the free energy takes the free-field 
form but with the one-particle energy cj(k) replaced by e(k), the latter satisfying a temperature- 
dependent integral equation with kernel related to a zero temperature form-factor of the trace 
of stress-energy tensor. For 2D integrable theories the approach reduces to the thermodynamic 
Bethe ansatz. For relativistic theories, a thermal c-function C qs (T) is defined for any D based on 
the coefficient of the black body radiation formula. Thermo dynamical constraints on it's flow are 
presented, showing that it can violate a "c-theorem" even in 2D. At a fixed point C qs is a function 
of thermal gap parameters which generalizes Roger's dilogarithm to higher dimensions. This points 
to a strategy for classifying rational theories based on "polylogarithmic ladders" in mathematics, 
and many examples are worked out. Other applications are discussed, including the free energy of 
anyons in 2D and 3D, phase transitions with a chemical potential, and the equation of state for 
cosmological dark energy. For the latter we show that — 1 < w < —1/3. 
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I. INTRODUCTION 



Quantum field theory at finite temperature and density was developed long ago mainly 
to deal with condensed matter systems at low temperatures. The subject was later taken up 
in the relativistic context where the interest was mainly in high temperatures Q, 3], with the 
aim for example to study hot QCD 4J. In both of these contexts the main computational 
tool is the beautiful but sometimes awkward sums over Matsubara frequencies in pertur- 
bation theory. It was later realized that to capture some important physics, such as the 
gluon relaxation rate, a re-summation scheme was necessary, beyond the re-summation of 
"ring diagrams", at least to re-sum the hard thermal loops[3j]. In spite of this important 
progress, finite temperature field theory remains difficult, especially when dealing with phase 
transitions, and it is worthwhile to explore alternative methods. In the condensed matter 
context other kinds of self-consistent re-summations are often used, for example in the finite 
temperature gap equation of the BCS theory of superconductivity. The main purpose of 
this paper is to develop a new approach to finite temperature field theory in any dimension 
D in a similar but essentially different spirit. Because the object we are dealing with, i.e. a 
quantum field functional integral in finite size geometry, is so general, there are many possi- 
ble applications; we will touch on a number of them, such as central charges at fixed points, 
classical statistical mechanics of Ising-like models in any dimension, anyons and cosmology. 

In sections II and III we introduce the usual quantities of quantum statistical mechanics, 
the free energy, pressure, entropy, etc. and express them as one-point functions of the stress- 
energy tensor. It is meaningful to express everything in terms of a temperature dependent 
c-function C qs (T) which at a renormalization group (RG) fixed point becomes the coefficient 
in the black body formula. (The subscript qs stands for quantum statistical and is meant 
emphasize the distinction with other c-functions.) We show that in 4D, this thermal c- 
function is different from other c-functions based on conformal anomalies [lj. We also argue 
that thermodynamic principles, such as positivity of the entropy, do not lead to any "c- 
theorem's" for C qs . 

C qs (T) is expected to be important for equations of state in the thermal universe. In 
section IV we relate the dark energy parameter w to a temperature derivative of logC qs . 
Using the entropy bounds on C qs we argue that for negative pressure — 1 < w < —1/3. 

The main workings of this paper begin in section V. The idea is to assume that the 
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partition function can be approximated by the standard free field form, but with the one 
particle energies replaced by a quasi-particle energy e(k). A saddle point approximation 
to the partition function then leads to an integral gap equation satisfied by e. The kernel 
in the integral equation is then related to zero temperature two-particle form factor of the 
trace of the stress-ener gy t ensor Tjf. This approach is modeled after the thermodynamic 
Bethe ansatz (TBA) 



32, 



34, 



, which is exact for the 2D integrable quantum field theories. 
For the non- integrable theories in D dimensions we don't expect the results to be exact, but 
rather to be a potentially useful approximation. In section VI we check that to lowest order 
our formalism gives the same result as the two-loop Matsubara formalism for the scalar 4 
theory. 

At an RG fixed point, C qs becomes a constant c qs which is an important characteristic of 
the conformal field theory. There are no known methods for computing c qs for interacting 
fixed points in dimension D > 2. As a step in this direction, in section VII we describe how 
to compute c qs in our formalism. The basic mechanism, which we believe to be general, is 
that at a fixed point the quasi-particle energy e develops a gap A that depends only on the 
temperature A = gT. The gap parameter g can be computed from an algebraic limit of the 
gap equation, and c qs is a function only of g which is a higher dimensional generalization of 
Roger's dilogarithm. In this way a fixed point theory is deconstructed into the thermal gap 
data, which depends only on interaction parameters that are an integral of the form factor 
of T». 

The classification of theories in any dimension with rational c qs is shown to be very 
closely related to the construction of polylogarithmic ladders in mathematics [3], in partic- 
ular classical analysis and number theory. Various aspects of this are explored in sections 
IX- XI. Given a rational c qs that can be obtained from the higher dimensional polylogarith- 
mic expressions of section VII, one can work back- wards by reconstructing the algebraic gap 
equation. The latter gives some information on the interaction parameters. For integrable 
theories in 2D, due to the tight constraints of integrability, one can in many cases entirely 
reconstruct the S-matrix of the theory. In higher dimensions we do not as yet have the tools 
necessary to complete the reconstruction of the fixed point theory. However we believe our 
analysis points to the existence of interesting conformal theories in higher dimensions and we 
can at least reconstruct some of their properties, such as particle spectrum and interaction 
parameters. 
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An interesting example worked out in detail is that of a single bosonic particle with gap 
parameter g = log 2. The polylogarithmic identities leading to c qs = 1/2 in 2D and c qs = 7/8 
in 3D where known to Euler and Landen. Though c qs = 1/2 is the central charge of a free 
fermion in 2D, the physical meaning of our bosonic theory and its possible connection to 
the Ising model will remain unanswered. The reason it has not been studied before is that 
previously nobody has been able to make sense of a bosonic quantum statistical mechanics 
in the TBA framework j^]. 

In section XII it is shown how a purely imaginary chemical potential can describe frac- 
tional statistics particles. In this construction it is necessary to work with gap parameters 
in the complex plane z = e~ 9 . For free "anyons" the gap parameters z are on the unit circle. 
In 2D the free anyons correspond to (3 — 7 systems in a continuous sub-regime of the spin 
a = d /2ti of the particles where < 1? < 2n. When the spin is rational, so is c qs . Extending 
this construction to 3D we propose a formula for the free energy of free anyons. Here c qs is 
only rational for special values of the statistics parameter such as djln = 1/3. 

In section XIII we reintroduce the chemical potential and analyze the phase structure. 
This leads to general statements concerning the existence of fixed points depending on the 
interaction strengths and chemical potential. 

II. FREE ENERGY AND THE STRESS-ENERGY TENSOR 

In this section we introduce the main physical quantities that are the subject of this 
paper; we mainly collect some well-known formulas. From the most general point of view, 
we will be dealing with a euclidean functional integral of fields in D = d + 1 dimensions 
where one of the dimensions is compactified into a circle of circumference (3. This object 
arises in a variety of physical contexts. It is the finite temperature T partition function of a 
quantum many-body system or quantum field theory in d spacial dimensions where the extra 
dimension is the euclidean time and /3 — 1/T. It can also represent a partition function 
of classical statistical mechanics in D spacial dimensions with one finite size direction of 
length (3. We will mainly be using the language of quantum statistical mechanics. Most 
formulas will be expressed in terms of (3, but some formulas will be expressed in terms of 
the temperature T when this makes the physical context more transparent. Some formulas 
will be expressed in terms of D, others in terms of d, depending on what is more convenient. 
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The partition function is 

Z(p,n) = Tre -^(H-^N) (1) 

where H is the many-body hamiltonian, N the particle number operator, and /j, the chemical 
potential. The free energy density, T = F/V, is defined as: 

where V is the <i-dimensional spacial volume. The energy density S per volume and pressure 
p are defined as 

l dlogZ g(OT \ d\ogZ 

V dp dp 1 V p dV [ ' 

The above quantities are related to the stress-energy tensor of the theory T MJ ,(x), /i, v — 

— * 

0, 1, .., d. The hamiltonian H and momentum operators P are given by 

H = J rfx T 00 (x), Pi = J dxToi(x) (4) 

where <ix = c? d x. The energy density and pressure are related to finite temperature one-point 
functions of T^ v : 

(Too)p = £, (Tu}(3 = p (5) 

(no sum on i), where denotes the finite temperature correlation function. (The one-point 
functions are independent of x by translation invariance.) The trace of the stress-energy 
tensor, = X^=o-^> can tnus ^ e elated to the free energy: 

{Tp p = {Pd p + d+l)T (6) 

In the sequel we will also deal with the entropy density S and number density n: 

The first equation is just the first law of the thermodynamics. 



III. QUANTUM STATISTICAL C-FUNCTIONS IN D DIMENSIONS 

For a free gas of bosons or fermions, the partition function is well-known in the limit of 
infinite volume V: 

/dk 
_ log (1 =F e-**) (8) 



where the upper (lower) sign corresponds to bosons (fermions), c?k = d d k, and is the 
single particle energy as a function of the momentum k = k, which need not be relativistic. 
For a massless relativistic theory with = V\s?, the integrals for the basic thermodynamic 
quantities can be done analytically. For example in 4D, the energy density is the well-known 
black body formula: 

£ = c — T 4 (9) 
q 30 v ; 

where c qs = 1 for a free boson and 7/8 for a free fermion. One can also relate the entropy 
S and the number of particles N = nV: 

S = e T N (10) 
where e_ (e+) is the entropy per particle for bosons (fermions): 

and ((z) is the Riemann zeta function. Finally the ideal gas law reads: 

pV = T+T NT < 12 > 

One question that will be addressed in this paper is how interactions can change the 
coefficient c qs in eq. JHJ) even at a fixed point. This leads us to introduce a quantum 
statistical c-function, C qs (T), which depends on temperature, and at a fixed point becomes 
in 4D for example c qs in eq. 0. This quantity can be defined in any dimension. If the theory 
is relativistic and conformally invariant, the trace of the stress tensor vanishes {T^)p = 0. 
Equation (JBJ then implies that T = const. / (3 d+1 . This leads to a natural definition of a 
thermal c-function based on quantum statistical mechanics: 

^-^T^,A) (13) 

where a d is a normalization constant. The function C qs depends on (3 and the couplings 
and masses of the theory, here denoted generically as A. Because the definition eq. (|13jh 



which is essentially motivated simr 

r 

appeared in a number of works j 



jy dimensional analysis, is so natural, it has previously 

the finite size, rather than finite 



imply Dy dimensional anaivs 



temperature, version was considered.) This paper is not so much concerned with proving 
or disproving c-theorems, but rather with developing a new framework for computing C qs . 
Nevertheless, some remarks on c-theorems can be found later in this section. 



The trace of the stress-energy tensor obeys renormalization group equations (RG) since 
it is a one-point function of a quantum field theory on a finite size geometry, the length scale 
being (3. Thus: 

Q(e s p, X(l)) = e- Ds Q(P, X(l + s)), Q((3, A) = (Tfi p (14) 

Above, / is the log of the length scale, and the / dependence of A is determined by the RG 
beta-functions: d\/dl = A(A). The above equation shows that C qs (/3) as a function of (3 
tracks the RG flow. Using eq. (0) one has: 

Q((3) =S-dp = a d T d ^ ^ (15) 

Thus at an RG fixed point G = 0, C qs is a constant independent of (3, and is an important 
characteristic of the fixed point theory. We will use the lower case c qs to denote the (3 
independent fixed point value of C qs . 

The entropy and energy densities in terms of C qs are: 

S = ?±{(d+l)-/3dp)C 9 (16) 
£ = ^{d-pd p )C 9 (17) 



At a fixed point one has: 



£ = d ~^m c ^ S = {d+ 1)^ c qs (18) 



so that S/S = ^T. 

We chose ad so that c qs = 1 for a free, massless relativistic boson. As we will show below, 
this leads to: 

a D ^ = k- d ' 2 T{D/2)C{D) (19) 

where V is the gamma-function and as before ( the Riemann zeta-function. When D is 
even, ((D) is rational up to powers of ir: ((D) = 2 d ~ 1 7t d \Be>\/ D\, where Bp are rational 
Bernoulli numbers. This is not the case for D odd; in fact it has only been proven relatively 
recently that ((3) is irrational. We checked numerically that ((D) /ir D for D odd doesn't 
appear to be rational. As we explain below, this appears to be related to the fact that there 
is no conformal anomaly for D odd. For D = 2, 3, 4 the expression for the free energy then 
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becomes: 

7TT 2 

= —Cqs(T) C -^- (D = 3) (20) 

7T 2 T 4 

= -C QS (T)- 

qsV ) go 

Since in recent years there has been renewed interest i n p ossible generalizations of 



(D 


= 2) 


(D 


= 3) 


(D 


= 4) 



n 



Zamolodchikov's c-theorem |l3j] to higher dimensions [rsl. Q, |20j, |2l|, we now discuss the 
relation of c qs to other proposals in higher dimensions. Let us first review the success story 
in D = 2. In a conformally invariant theory, c can be defined by the leading singularity 
of the operator product expansion of two stress-energy tensors, or equivalently the central 
extension in the Virasoro algebra[14]. The same central charge c appears as a conformal 
anomaly in the conformal field theory: 

(TP = -£* (21) 

where R is the scalar gravitational curvature. Using the above relation one can show that c 
governs the finite size effects 3, till of the partition function and thus at a fixed point c is 
the same as c qs . 

For the purpose of studying the RG flow between fixed points there are two known 
alternative approaches to defining a scale dependent c that tracks the flow and equals the 
Virasoro central charge at the fixed point. One is to define C qs (/3) as abovej^]. The other is 
the content of the c-theorem where one defines a function c z (L) from the two-point functions 
of Tjjj, at zero temperature but at finite separation L. The c-theorem states that c z (L) 
monotonically decreases as a function of the increasing length scale L. There is sometimes 
the misconception that C qs and c z are essentially the same. This misconception has led 
to the expectation that C qs also always decreases, even in 2D. However these two c's are 
potentially very different, as one is related to a one-point function at finite temperature, 
the other a two-point function at zero temperature. The thermodynamic foundation for C qs 
actually makes it a much richer physical quantity, and obviously one that can be defined in 
any dimension. 

Though some discussions of the c-theorem vaguely invoke notions in statistical physics 
like irreversibility and loss of information, in reality c z has no statistical mechanical meaning. 
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The intuitive understanding is more accurately that c z is a measure of the massless states 
of the theory, and since massive states decouple at low energies, c z decreases [l^. Ironically, 
whereas C qs (f3) does have a thermodynamic meaning, no theorem has ever been proven for 
it. Indeed, though C qs is known to decrease for theories with UV and IR fixed points, it can 
actually oscillate in physically sound theories with RG limit cycles 0, 24 1. 

As eq. (115)1 obviously shows, a c-theorem for C qs follows from positivity of G: if G > 0, 
then dC qs /dT > 0, i.e. C qs increases with temperature. Let us refer to the above condition 



as Q-positivity. In 11] many interesting examples were presented that exhibit the expected 
property c qs IR < c qs uv ; however examples where C qs does not decrease monotonically with 
decreasing temperature were also discussed. 
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24j . Suppose 



A general manner in which Q-positivity can be violated is as follows 
a quantum field theory is described by a conformally invariant action S'cft perturbed by 
operators O a \ 

S = 5 cft + ! ^ X aO a (22) 

where A^ are positive couplings. Let \a be the beta function for A^ and Ta the scaling 
dimension (including anomalous) of O a . Then the beta functions A^ = d\A/d\ogL, where 
L is a length scale, to lowest order are 

X A = (D-T a )Xa + 0(X 2 ) (23) 

The main point is that the trace of the stress-energy tensor receives quantum corrections, 
and since it must be zero at a fixed point where the beta functions are zero, one must have: 

A 

The above formula is well known and is easily verified in 2D to lowest order in conformal 

n 

perturbation theory [25]. Consider first a theory that can be formulated as a perturbation 
of an ultra-violet fixed point by relevant operators, which implies Ta < D. Then the beta- 
functions are positive to lowest order. Furthermore, for relevant perturbations, because of 
the anomalous dimensions of the couplings, there is often no higher order corrections to 
the beta functions since higher powers of A^ do not have the right dimension. So in this 
situation, is generally positive and the c-theorem holds. 

The above arguments clearly point to the way in which the c-theorem can be violated. 
If the O a are marginal, Ta = D, and the beta functions start at 0(A 2 ), and there are no 
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general constraints on the sign of the beta function. This was suggested as the origin of the 



23, 



2J , though the 



violation of the c-theorem for the 2D models with limit-cycles in the RG 
issue is not entirely resolved. 

What constraints on the flow of C qs follow from the most basic thermodynamic principles? 
Note that the entropy itself, and not its change, already depends on the flow of C qs , eq. 
(HHJ). One constraint is that the entropy density should be positive. Since the pressure is 
proportional to C qs , there are two cases depending on the sign of C qs . Using eq. (fTfij) and 
(fT7jl (in terms of temperature T): 



positive pressure 



S > 



£ > 



9 log C qs 

«91ogT 

9 log C qs 

<91ogT 



> -d- 1 



> -d 



(25) 



negative pressure : S > =>- - ^ ^ qs < —d — 1 (26) 

a log T 

c n <9 log Cqs 

o log J 

Note that in the positive pressure case, a positive entropy could have a positive or nega- 
tive energy density. On the other hand in the negative pressure case, a positive entropy 
necessarily implies a positive energy density. 

Let us now turn to the conformal anomalies in higher dimensions. In D dimensions 
has scaling dimension D and the curvature R has dimension 2. This already indicates that 
there is no conformal anomaly for D odd since no power of the curvature has dimension D. 
As mentioned above, this appears to be reflected in the formulas eq. (J2*U|) by the irrationality 
of ({D) up to geometrical powers of 7r for D odd. Again the quantity C qs presents itself as 
the most natural analog to c for all even and odd dimensions. 

For D = 4 the conformal anomaly reads 0, |3| 



where W is the Weyl tensor and R the dual curvature. (We have rescaled c, a so that 
c = 1 for a free boson.) Cardy conjectured an a-theorem for the anomaly a, and recently 
there have been many interesting examples of supersymmetric gauge theory supporting 
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the conjecture [it], |20|, |2l|. Many of these examples rely on Seiberg duality. It should be 
emphasized that the above equation is valid only at a fixed point, and thus does not provide 
a definition of scale-dependent c, a that interpolate between fixed points. Since c and a are 
related to one-point functions of Tff, they are analogous to c qs rather than c z . However 
one can easily show that no linear combination of c,a equals c qs . Consider a free theory 
consisting of Nq scalars, spin 1/2 Majorana fermions, and N\ spin 1 gauge bosons. The 
anomalies are known [2l|, |22|: 

c = N + 3N 1/2 + UN! 

a = (2N + llN 1/2 + 124iVi)/6 (28) 

On the other hand c qs only depends on the statistics of the particles and their number of 
degrees of freedom, which are 2 helicities for the fermions and 2 polarizations for the gauge 
bosons. Using c qs = 1 for a boson and c qs = 7/8 for a fermion: 

c qs = N + 7 -N 1/2 + 2N X (29) 

Noting that 3(c~— a)/2 = N + 7A^!/ 2 /4 — 13iVi, one concludes that in a free theory of only 
scalars and spin 1/2 particles, c qs = 3(c — a)/2, however for higher spin particles there is no 
linear combination of c~and a that equals c qs . 

Since the equation (|2*Tj) involves the curvature, the anomalies a can determine some 
finite size effects. However unlike the D = 2 situation, a apparently do not govern the 
finite size effects having to do with finite temperature. Indeed, lattice simulations of finite 
temperature pure QCD observe the behavior eq. © at high temperature with c qs = 16 
which is consistent with an SU(3) gauge theory with 8 gauge bosons The fact that c qs 
is here given by the free field values is an indication of asymptotic freedom of QCD [J]. 



IV. QUANTUM STATISTICAL C-FUNCTION IN COSMOLOGY 

The free field fixed point values of C qs are used extensively in studies of the thermody- 
namics of the hot early universe. For instance the energy density is taken to be given by the 
black body formula © for a radiation dominated universe. This is clearly an approximation 
in a number of respects. First of all, even at a fixed point, the interactions can change the 
value of c qs , unless the theory is asymptotically free. This will be the subject of section 
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VII of this paper. Secondly, the universe is probably not at a fixed point. Our comments 
on cosmology will be confined to this section, where we primarily phrase some cosmological 
properties in terms of C qs and use this to place bounds dark energy that is not a cosmological 
constant. In other words, we place bounds on dark energy that is hypothetically made up 
of particles. 

There is growing observational evidence for an acceleration of the expansion of the uni- 
verse. The cause of this phenomenon is commonly referred to as dark energy. (See for 



instance, |27 



eter w: 



28, 



291].) Dark energy is usually characterized by the equation of state param- 



One of the Friedmann equations expresses the acceleration in terms of S,p: 

-(S + dp) = — [- + d)p (31) 



R 3 y w 3 \w 

where throughout this section d = 3 and R is the cosmic scale factor. Accelerated expansion, 
i.e. dark energy, corresponds to R/R > 0. If the pressure is positive, then dark energy 
requires — 1/d < w < which corresponds to a negative energy density. If on the other 
hand the pressure is negative, dark energy requires w < —1/d with positive energy density, 
or w > but with negative energy density. 

The parameter w can be expressed in terms of C qs : 

1 , <91ogC as 

- = d+ S ^ s 32 
w a log T 

At a fixed point, w = 1/d regardless of the value of c qs . However away from a fixed point, it 
can be different. Already one sees that if dark energy is matter, it is currently not at an RG 
fixed point. The thermodynamic bounds coming from the requirement of positive entropy, 



eqns. f|25!26|) . place bounds on w: 

positive pressure : S > w < — 1 or w > 

negative pressure : S > — 1 < w < (33) 

Thus one sees that positive pressure dark energy, which required —1/d < w < 0, is incon- 
sistent with positive entropy. More importantly, taking into account the entropy bounds, 
negative pressure dark energy necessarily corresponds to— 1 < w < — 1/d. 
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The standard candidate for dark energy is a cosmological constant A, where the stress ten- 
sor T^ v oc Kg^ v . In Minkowski space, this can correspond to the stress-energy tensor of vac- 
uum energy. From the Lorentz structure of the flat space metric 77^ = diag(l, —1,-1,-1), 
this necessarily implies p = —£ and w — — 1. Interestingly, the lower limit w — — 1, which 
we obtained from thermodynamic principles, precisely mimics a cosmological constant. 

If dark energy is not a cosmological constant, then in terms of C qs the most likely scenario 
for dark energy is that the universe is currently at a temperature where C qs is negative with 
dCqs/dT > 0, i.e. C qs is decreasing as the temperature continues to decrease. The fate of 
the universe then depends on whether there is an RG fixed point at zero temperature. If 
there is a fixed point, then as explained above, w will eventually become positive and equal 
to 1/3. 

As we will see, obtaining a negative C qs is not so straightforward in our formalism. One 
can study this by examining fixed point values c qs analytically and see what is required for 
c qs to be negative, since C qs is bound to be negative in the vicinity of the fixed point. As we 
will show in section XIII, this seems to require an imaginary part to the chemical potential. 
One interpretation of an imaginary chemical potential in 2D and 3D will be given in section 
XII, i.e. as corresponding to fractional statistics particles. Analogous Hopf terms in 4D 
would require Im{p) = iirT since there are only bosons and fermions in 4D; one example is 
skyrniionsj^c]]. But there may be other interpretations to an imaginary part, for example a 
relaxation rate. 

V. QUASI-PARTICLE RE-SUMMATION 

This section is rather technical, so let us begin by drawing attention to the main results: 
eq. ()42)1 is the desired form of the partition function, where e satisfies the integral gap 
equation ()4()j) . The kernel appearing in this integral equation is given by the form factor eq. 

Interacting field theories, even at finite temperature, normally have a particle description. 
In many complex condensed matter systems one is accustomed to dealing with quasi-particle 
excitations of definite statistics, for example in Landau's theory of a Fermi liquid for metals. 
The free particle form of logZ, eq. (jHJ) is appealing since it incorporates the bosonic or 
fermionic statistics of the particles. The eq. (JHJ) can of course be obtained from a functional 
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integral in finite size which involves the Matsubara frequencies, but only with some effort; 
the simplicity of the final answer suggests there is a simpler, particle description. Based 
on these observations, we formulate a quasi-particle description of Z with the following 
properties. We suppose that the free particle form eq. (jHJ) continues to hold but with 
replaced by a quasi-particle energy e(k). All the effects of the interactions, and additional 
finite temperature effects are all contained in e. 

There are no general arguments showing that finite temperature field theory can be 
organized in this fashion. The formulation should be viewed as a potentially useful approxi- 
mation that can capture some physical effects in a transparent way in comparison with other 
approaches. However for 2D integrable theories this formulation is exact and known as the 
thermodynamic Bethe ansatz, and this certainly provided some motivation and insights. 

We first describe a particle description of the partition function. To proceed, introduce a 
density of particles p(k) so that p(k)<ik represents the number of particles with momentum 
between k and k + rfk. Similarly introduce the level density pi(k) which counts the available 
quantum states. In a free theory pi just counts the states for particles in a box of volume V 
and pi = V I (27i) d . In an interacting theory the densities p and pi are clearly not independent: 
when a new particle is added it affects the available levels because of the interactions. 
For integrable theories in 2D the quantization condition on allowed momenta leads to a 



relation between p and pi involving the 2-particle factorized S-matrix[34, 35]. (See section 
XB.) It is not clear that this kind of quantization condition can actually be formulated in 
higher dimensions. So let us for now simply assume that there is a linear relation between 
the densities that depends on the momenta. This can generally be written in terms of a 
temperature independent Fredholm integral operator (kernel) G: 

V f dk' 1 
"'( k) =(2^ + /(2^"( k ' )G ( k ' k ^ < 34 » 

The kernel G will be determined below. However as we now show, any choice of G leads 
to the desired form of the partition function. Retaining the chemical potential, one has 

log Z = -13 J dk (co k - /i)p(k) + $ (35) 

where $ is the entropy. The entropy depends on the statistics of the particles: 

$ = ± J dk((pi±p) log(pz ±p) - pi log pi =F p log p) (36) 
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where the upper/lower sign is for boson/fermion statistics. (In the subsequent formulas the 
boson case will always be the upper sign.) We now look for a saddle point approximation 
to Z. Varying with respect to both densities: 

51ogZ = {-(5{u -fi) + \og( Pl /p ± 1)) 5p ± log(l ± p/ Pl )5 Pl (37) 

The saddle point is subject to the constraint eq. relating p,pi. The saddle point 

equation 5 log Z = then has a simple form. Define the filling fractions 

/ = - (38) 
Pi 

and parameterize them as follows: 

h = — (39) 
Then the saddle point equation becomes a non-linear integral equation involving G: 

e(k) = (o; k - fi) ± i J ^ G(k, k') log (l =F e'^) (40) 

where we have defined: 

dk = , ^ (41) 

= (2vr)^ k 1 ; 

which is Lorentz invariant when the theory is relativistic. Plugging this back into log Z one 
finds that it has the desired form: 

logZ = TV Iwr log (1 T e ^ £(k)) (42) 

We will refer to the equation (J4*U|) as the integral gap equation. 

The fundamental thermodynamic quantities, obtained by differentiation, have simple 
expressions: 

r rlV 

£ = J ^y d f T (k) (1 + (3d,)e (43) 
dk „ s de 



n = -JWf^T, (44) 

Also, using Vk - k = d, one has for the pressure 

P=1 dlwr f ^ k)k - V ^ k) (45) 

Using the integral equation satisfied by e, one can obtain an integral series representation 
for the above quantities. Introduce the shorthand notation 

(G*0)(k) = y^G(k,k')0(k') (46) 
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for an arbitrary function 0(k). (Where there is no chance for confusion we will sometimes 
drop the *.) Then using the integral equation (|40|) one can show: 

(1-K)(l+Pdf,)e = u (47) 

where K is the kernel 

if(k,k') = G(k,k')/ T (k') (48) 

Similarly 

(l-K)9 Ai e = -l (49) 
Using now (1 — K) _1 = J2m=o K m one obtains 



m=0 



n = E/7^d/T(k)K»*l (51) 

m=0 J ^ ' 

So far, our equations are quite general and valid whether the underlying theory is Lorentz 
invariant or not. For the remainder of this article, we assume the theory is Lorentz invariant 
with: 



uj k = Vk 2 + m 2 (52) 

We can now derive an expression we will need for the pressure. Introduce the 4-vectors 
k^ = (u?k, k) with k^kn = m 2 . The kernel G can only be a function of the Mandelstam 
variable: 

s(k,k') = (k^ + k' t ,) 2 = 2(m 2 + k%) (53) 

(The other possible variable t = (k^ — k'^) 2 = 4m 2 — s is not independent.) We will need 

2 2 
Vk s = — (t^k'k — CG> k k'), Vk' s = — (u)ek- — ^k'k) (54) 

The above implies for the kernel G: 

u h V k G(k, k') = -uv V k ' G(k, k') (55) 

Using the integral equation for e one can now show 

(1-K) V k £ = V k cu (56) 
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where K (k, k') = uj^ 1 K(k, k')uv- This leads to 

dk 



P 



d ^ 

m=0 



(2tt)« 



/ T (k) k-K m * V 



10 



(57) 



What remains is to determine the kernel G. Consider the one-point function (Tff)p. From 
first principles one has (// = 0): 

1 



TO/ 



; Tr (T£e- /3H ) 



(58) 



The left hand side can be computed from the free energy eq. © whereas the right hand 
side can be computed by inserting a multi-particle resolution of the identity into the trace 
and expressing the result in terms of the form factors of T^. In this way we will relate the 
kernel to a certain zero temperature form factor. This computation was carried out for 2D 
integrable theories in [3^], and we will subsequently follow this reference closely. 

Let us first determine (T£) from the free energy. Since {Tff)p = £ + (D — 1) p, using eqs. 
(150157)) one obtains: 

dk 



m=0 J ^ > 



(59) 



Next we describe a form-factor computation of (Tj£)p. Introduce particle states |k) nor- 
malized as follows: 

(k'|k> = (2tt)V 5 {d \k-k') (60) 
Inserting a multi-particle resolution of the identity in eq. ()58|) one finds 



n=0 J 



8=1 



(k n ...k 1 |T^|k 1 ..-k n ) 



(61) 



The form factors (k n ■ • • |T^|ki ■ ■ •) are the zero temperature form factors. These form factors 
have two kinds of contributions: disconnected contributions proportional to 5 functions and 

nn 

the connected part which involves no 5 functions. As explained in [2JJ, |22j , the effect of the 
disconnected parts is two-fold: some cancel the overall 1/Z and others convert the e - ^ into 
the filling fractions / T . Because the partition function involves the quasi-particle energy e, 
the filling fractions / are also the quasi-particle energy ones given in eq. (|46j). The final 
result is then: 



TO/ 



co 1 

T- 

n=0 



i=l 



(k n ---k 1 |T^|k 1 ---k n ) c 



(62) 
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Finally comparing eq. ()62|) with the n-th term in the sum eq. (|59j) one finds 

(k^- • -ki|T^|k a • ■ -k n ) conn = (uj kl u kn - k x ■ k n )G(k n , k n _i) • • •G(k 3 ,k 2 )G f (k 2 ,k 1 ) +perm. 

(63) 

The above factorization is exact for the 2D integrable theories as a consequence of factoriz- 
ability of the S-matrix; for non-integrable theories it is likely to be an approximation. The 
above equation implies: 

(k|T;|k) conn = m 2 (64) 

which is a normalization condition. The kernel is then identified with a two-particle to 
two-particle form factor: 

G(k, k') = ( k) k| T M| k ' ) k)conn (65) 

where 2k^k' = s — 2m 2 . 

In the sequel we will need a scaling equation obeyed by the kernel. G depends on the 
momenta k and couplings and masses denoted as A with scaling dimension d\. Using only 
the fact that G has scaling dimension (1 — d) one has: 

G(e% e l k'; e^A) = e {l ~ d)l G(k, k'; A) (66) 

It is straightforward to extend the above formulas to many particles. Let there be N 
particle species, each with statistical parameter s a = ±1, where s a = 1 corresponds to a 
boson, and m a their masses, a = 1, .., N. The equations become: 

log Z = -V J j^r d s a log (1 - s a e~^) (67) 

a 

e a = K, a - 1^) + - J2 G ab * s b log (1 - s b e-? £b ) (68) 
" b 

where u^ a depends on m a . The kernels are 

G a ,(k,k') = - J— (k,a;k',6|T;|k',6;k,a) (69) 

lK a K vf> 

In summary, the quasi-particle approach we have developed is implemented as follows. 
One first computes the kernel G from the zero temperature form factor eq. ()65j1 to whatever 
order in perturbation theory is feasible. One can then solve for the quasi-particle energy e 
by numerically solving the integral equation (|40|). Finally the partition function is given by 
the simple eq. (|4~2"|). 
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VI. COMPARISON WITH PERTURBATION THEORY IN THE MATSUBARA 
FORMALISM. 



The standard methodology of finite temperature field theory involves perturbation theory 
where one must sum over Matsubara frequencies. Since this looks completely different from 
the approach of the last section, it is important to check that both approaches agree at 
least to lowest order. In this section we perform this check for the 4 scalar field theory 
in any dimension. The lowest order contribution to the free energy, eq. ()86|) below, is 
well-known In the latter formalism, this lowest order correction is a 2-loop vacuum 

Feynman diagram. As we will see, to obtain this standard Matsubara result in our formalism 
is somewhat subtle. 

The interacting scalar theory is defined by the lagrangian 

£ = - ^ 2 - ^ (70) 

The usual stress-energy tensor is 

T ia , = d ft <f>d v <l>-ri lu ,£ (71) 

The problem with the usual stress-energy tensor, is that even in the conformal limit where 
the theory is free and massless (A, m = 0), the above stress-energy tensor is not traceless. 
It is known in this limit how the stress-tensor can be improved so that it is still conserved 
but now traceless 

T- w = Tfu/ + (v,<9 a (<R0) - <9,(<R0)) (72) 



In the interacting case, the further improvement eq. (|24|) is necessary: 

TjH = ^ (73) 

where A is the beta function for flow toward the infra-red. (A differs by a sign from the usual 
convention of flow toward the ultra-violet in high energy physics.) Since the field has mass 
dimension (D — 2)/2 in D dimensions and the operator </> 4 has dimension 2(D — 2), 

A = (4 - D)X + 0(A 2 ) (74) 

We can now compute the form factor that determines the kernel. Expand the field in the 
usual way: 

ct>U) ~- 



a^e + a^e I (75) 



(27r) d v/2^ 
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where 

[a k/ ,al] = (2^ d 5^(k-k') (76) 
The states |k) with the normalization in section IX are then: 

|k> = 4|0) (77) 

Expanding the fields <fr i n eq. (J75|) and evaluating the matrix element one finds: 

(k,k'|T;|k',k) = ^ (78) 

Thus the kernel G to lowest order is 

G < k ' k '> = 4(^) (") 

where s(k, k') is the Mandelstam variable eq. (|53p. 

The integral gap equation implies that the lowest order correction to e is: 

e « uj u + * log(l - e"^) (80) 

Inserting this into the free energy, one finds 

jr = jr + jr 1 + o(A 2 ) (81) 

where 

= ^ y ^1 ^ /o(k) u; k G(k, k') log(l - e-^') (82) 
Above, the filling fraction is the free-field bosonic one: 

= ( 83 ) 

Next we express the kernel as a derivative. For simplicity we only consider the massless 
case. One can easily show: 

^k.V t , fi)=i, -k-.Vfc fi)=i (84) 
V s / s V s / s 

Using these two equations in turn, and integrating by parts, one has: 

^ = ~\f §kj §L /oOO/oOO ^ (85) 



{D-2)T 1 + jj^f^ /o(k)/ (k') 
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The extra (D — 2)T\ comes from Vk • (k/fc) = (D — 2)/k. Adding the two above equations 
one obtains 

•F. = 5 f / «(k)) (86) 

nn 

The above equation is the well-known lowest order correction to the free energy Ul,|2|. Note 
that in 4D one needs to keep the leading term (4 — D)X in the beta function (J74"j) through 
to the last step in order to obtain eq. (jBEJ). 



VII. THERMAL GAPS AND c qs FOR INTERACTING FIXED POINTS 

In this section we analytically study the integral gap equation for relativistic theories 
near an RG fixed point. In this way we will obtain expressions for c qs for an interacting 
fixed point. Throughout this section the chemical potential /i = 0. For simplicity we first 
treat a theory with a single particle. 

In the approach developed in section V, C qs is given by the formula: 

fid r ju 

^^/SF^^ (87) 

As we argued in section III, at a fixed point C qs is a constant independent of (3. At a fixed 
point the theory is expected to be massless. This can occur in the ultra-violet (UV), i.e. at 
high temperature or high energy compared to the mass where ~ vk^ = k. We will focus 
on such UV fixed points, but some of our analysis can also apply to infra-red fixed points 
at low temperature since such a theory is expected to be massless from the beginning for all 
k and again = k. The fixed point properties are then determined by the solution to the 
massless integral gap equation: 

s = k ± —G * log (l =F e _/3e ) (88) 
P 

At a fixed point, G should depend only on dimensionless parameters A with d\ = 0, so that 
eq. (jHHjl implies 

G(e%e l k') =e ( - 1 - d)l G(k,k') (89) 
Then from eq. (jHHj) one can show 

e (e'k, 0) = e l e(k, e l (3) (at a fixed point) (90) 



21 



The fixed point value of C qs is then: 

CM = t— ! 4k k log (1 T e~ P£ ) (91) 

&d J 

where now e satisfies the massless integral gap equation (J88|) . By rescaling k in the above 
equation one shows that C qs (e l P) = C qs (/3). In summary, if e satisfies the massless integral 
gap equation (|SS]l. and G the scaling equation (|85j) then C qs in eq. (|57|) is a constant c qs 
independent of (3. 

An important feature of the massless integral gap equation (J88)) is that e generally de- 
velops a thermal gap A oc T at high temperatures. When k is small compared to the 
temperature T, the fc-term on the RHS of (jHHj) can be dropped and (3e is approximately 
constant. Let us define the dimensionless thermal gap parameter g so that: 

gT = A, fc<T, T large (92) 

For much higher than T then e ~ k. These features are shown in figure 1. In the limit of 
T — >• oo, e = A for all k, and g is a solution to the algebraic gap equation: 

g = ±h\og(lTe- 9 ) (93) 

where h is the dimensionless parameter: 

f dkG(0,k) (94) 



Henceforth, by gap equation we will mean the algebraic equation (}9"3*j) rather than the integral 
gap equation (|40jl. The interactions are encoded in the parameter h which we will refer to 
as interaction parameters. 

The quantity c qs should only depend on the gap parameter g. It can be computed to a 
very good approximation, if not exactly, as follows. In eq. ()87|) we wish to trade the integral 
over k for an integral over e. Differentiating eq. (jHHjl with respect to k one finds 



k = kd k e =F ^kd k G * L (95) 

P 



where we have defined L = log(l =F e /3e ). Inserting this into eq. (JHZJ): 



c qs = =F^- J (kd k e =f ^kd k G * L^j L(k) 



(96) 
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FIG. 1: The thermal gap. 
Using the property eq. (J55)l and a few integrations by parts one obtains: 



r (iS x / " de k d ~ l (±L(e) + ' 



1 a \ e^Tl 

When e > A, k ~ e, so we approximate k by e in the above equation. The final result is: 

1 



(97) 



c qs (g) 



rp-i)cp) 



Ms) 



(98) 



where 



3 (ITT) f d * l0E(1 * e ~> + 

(We have used f dk = Qd f k d ~ l dk, Vt^ = (2n d ^ 2 /T(d/2)), and the duplication formula for 
the T-function.) 

The above integral can be expressed in terms of the the polylogarithm, or Jonquiere, 
function: 

00 

Li «W = £i (100) 



k=l 



For D = 2, 3, 4 one finds: 



C 2 (g) = ±|Li 1 (±e^)±Li 2 (±e- 9 ) 

C 3 (g) = ±£lA 1 (±e-*)±gLi 2 (±e-°)±Ua(±e-s) 

U{g) = ±^U 1 {±e~ 9 ) ± g 2 U 2 (±e~ 9 ) ± 2gU 3 (± e - 9 ) ± 2Li 4 (±e- 9 ) 



(101) 
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Note that Lii(z) = — log(l — z). 

Since the argument of the polylogarithms is always ±e~ 9 , it is convenient to define the 
gap variable: 

z = e' 9 (102) 
To simplify further the subsequent notation, introduce the functions: 

r\ ~ nl 



Lr n (z) = U n (z) + V tJL\og r \z\U n „ r (z) + (-)" il?_l2 log- 1 \z\ log(l - z) (103) 



r=l 

Comparing with eq. (|101|) . one finds 



C d+1 (g) = ±T(d) Lr d+1 (±e- 9 ), (104) 



One then has the simple formula for real z: 



c qs( 2 ) = 77J-, t\ sLr d+1 (sz) = c D (z,s) (105) 
C(«+ 1) 



where as before, s = +/— corresponds to bosons/fermions. 
For many particles the gap equation becomes 



where 



In terms of the z a , 



b 

The quantity c qs then becomes a sum: 



9a = Y, h <* s 6 log(l-s 6 e- 96 ) (106) 

b 

h ab = J dkG ab (0,k) (107) 

l[(l-s b z b )- KbSb (108) 



a 



(109) 



Let us summarize the results of this section. The field theory has been deconstructed into 
the interaction parameters h ab , which determine the thermal gap parameters z a through the 
gap equation. The fixed point value of the central charge c qs is then a function of the z's. 

Since h ab determine the fixed point properties, they must be an RG invariant functions 
of the couplings. For the 2D integrable theories, it is well known how the fo's arise from 
coupling constants in the exact S-matrix. (See section XB.) Since we have relaxed the 
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integrability, let us illustrate how to compute h to lowest order in the 4 theory in 2D, 
which is not integrable. The kernel G was determined to lowest order in section VI. The 
interaction parameter h to lowest order is then: 

A f dk 1 . , 

h = / 7 110 

The above integral is convergent, and using A = 2 A + to lowest order one finds 



Since A has dimension 2, the ratio A/m 2 is a dimensionless coupling constant, and thus so 
is h. The 4 model may be considered as the weak coupling limit (b — > 0) limit of the 
sinh-Gordon model defined in eq. ()133J) where h oc b 2 . However the integrable sinh-Gordon 
model actually has h = 1 indepependent of b (see below) which shows the importance of the 
higher order corrections. 

For higher dimensions it is not clear how the structure of the kernel G can lead to finite 
RG invariant /i's, and this will remain an important unanswered question in this paper. 



VIII. CLASSICAL STATISTICAL MECHANICS: ISING-LIKE MODELS 

In this section we consider the application of the above formalism to classical statistical 
mechanics in D spacial dimensions. Since the physical context here is very different, the 
subscript c i will be used to denote quantities in classical statistical mechanics to avoid the 
likely confusion. The meaning of c qs in this context is the following. Consider a statistical 
mechanical system on a D dimensional lattice where all dimensions of the lattice are infinite 
except for one of length L, i.e. the volume of the lattice is VL where V is a D — 1-dimensional 
volume and V —>■ oo. At the critical classical temperature T c , the partition function should 
only depend on L. If periodic boundary conditions are imposed in the finite size direction, 
then L is the same as a finite inverse temperature f3. Thus the results of section III imply: 

lim ^-logZ.^c^ao^L 1 - (112) 

Consider an Ising-like model onaD dimensional square lattice. At each site i there is a 
spin degree of freedom Oi = ±1. The partition function and classical energy E are 

config. <i,j> 
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where < i, j > denotes nearest neighbors, and f3 c \ = 1/T C \. For D = 2 the usual Ising model 
is equivalent to a relativistic quantum field theory of a free Majorana fermion in euclidean 
2D. Though this is well-known, it is very non-trivial and relies on mappings to the XY spin 
chain and a Jordan- Wigner transformation. 

No higher dimensional version of this is known. We can explore the possibility of a 
higher dimensional version using the gap equation. In a disordered phase, since the spins 
are disordered, one expects that < E >= 0. From the first law of thermodynamics, one then 
obtains that the free energy is proportional to the entropy: F d = —S c \/f3 c \. Now suppose 
that at the critical point log Z c \ is equivalent to a field theory functional integral at an ultra- 
violet fixed point. Let us clarify that the temperature of the field theory T and the classical 
temperature T c \ have nothing to do with each other. As is well-known in the 2D case, the 
classical temperature appears as a coupling constant in the zero temperature field theory, 
the mass of the fermion being proportional to T c \ — T c where T c is the critical temperature. 
Using the results of section VII, one finds in the ultra-violet limit of the field theory for the 
entropy: 



where g is the thermal gap. Since V/(2ir) d is approximately a level density, then S ~ 
± l°g(l T e~ 9 ). From this equation, one sees that g = log 2 in the bosonic case is special 
since then the entropy resembles S = ^ log 2, which suggests a model with 2 states per site, 
such as an Ising model. 

Though the above discussion is very rough, g = log 2 indeed leads to the same central 
charge c qs = 1/2 as for the 2D Ising model (|98jl. For 3D, the same formulas again give a 
rational result, c qs = 7/8. This rationality is quite remarkable given the complexity of the 
various functions involved. The existence of rational theories will be explored in the next 
section. In 4D, rationality ceases and one finds c qs = Lr 4 (l/2)/£(4) = 0.97805427.... Clearly 
more work is needed in the 3D case, for example one could check whether for an infinite slab 
of thickness L, \ogZ c \/V scales as 7C(3)/167rL 2 . 

IX. RATIONAL THEORIES 

In 2D, theories with a rational c qs are of particular importance. For instance unitary 
theories with < c < 1 have been classified and are the so-called minimal models with 




(114) 
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rational c = 1 — 6/(n + 2)(n + 3), n = 1,2,... 14, 39]. Also, all known integrable scattering 



theories lead to thermodynamic Bethe ansatz equations with a rational c qs in the ultra-violet 
limit. 

The function Lr 2 (z) is the Roger's dilogarithm. Remarkably, the higher order functions 
Lr„ have also appeared in the mathematical literature [43]. They were not originally defined 
in terms of the integral of section VII. Rather, these functions arose as the unique func- 
tions that remove logarithmic terms in certain functional equations satisfied by Li n (z). We 
can now formulate the following problem. A rational theory in D dimensions requires the 
following identity: 

N 

J2sa L *D(s a Z a )=1l((D) (115) 

0=1 

where s a — ±1 corresponds to the statistics and z a are numbers that characterize the thermal 
gaps g a , z a = e~ 9a . If 1Z is rational, then c qs = 1Z is rational. We will refer to a relation 
of the form eq. (|115|) as a Rational Statistical Polylogarithmic Evaluation (RSPE). Every 
RSPE may be viewed as a non-trivial result in number theory. 

Of course, a relation of the form eq. (jll5|) does not guarantee the existence of a field theory 
with this c qs in the UV. However one can try to go as far as possible in reconstructing some 
basic data of a possible theory. Given the z a from an RSPE, one can use the gap equation to 
determine the interaction parameters h a ^. There then remains the difficult problem of finding 
a consistent field theory leading to kernels G a b that yield the h^. In 2D, if an integrable 
theory exists with c qs , the theory can often be entirely reconstructed, i.e all the scattering 
amplitudes can be deduced, due to the tight constraints of integrability. The encouraging 
fact is that, as we show below, we have been able to reconstruct the field theories for all 
known one-particle (N=l) RSPE's in 2D. We have also been able to reconstruct a 2-particle 
theory in 2D. 

The classification of fixed points of field theories is then closely related to the classification 
of RSPE's. An RSPE is a subclass of the more general polylogarithmic ladders constructed 
in 43j] . The latter do not have the restriction that the statistical parameters s a enter the 
relation as in (jll5|) . though some polylogarithmic ladders can be brought into the form 
of RSPE's using certain functional transformations (see below). The gap equation ()108|) 
is essentially what is called the cyclotomic equation in the construction of polylogarithmic 
ladders. The forerunners of dilogarithmic ladders (d=l) were first found by Euler and 
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Landen in the 1700's [4lL 1421 ] . New dilogarithmic ladders were not found until the 1930's by 
Coxeter and Watson. In the 1980's it was realized by Lewin that there is a vast extension of 
these previous few results. Despite this progress, polylogarithmic ladders are still considered 
somewhat mysterious and there are no classification theorems. Some of them have only 
been verified to high precision numerically. Ladders have even been discovered to order 
D = 17 



4611 . In the sequel we will present a number of RSPE's that were deduced from 
results in [43, 44, 45^ although they do not all appear there in the form they are given here; 
the diligent reader can easily verify them numerically (as we did). 

Let us first describe the only results known that are valid for all dimension D. Polyloga- 
rithmic ladders are constructed using functional relations satisfied by Lr n (z), without doing 
any explicit integrals. For any n, Lr n (l) = ((n). This corresponds to a free boson in any 
dimension with c qs = 1. The functions Lr„(z) satisfy a number of functional equations that 
are valid for any n. The duplication relation reads: 

Lr n (z) + Lr n (-z) = 2 l - n hi n {z 2 ) (116) 

From this it follows that — Lrd+i(— 1) = (1 — 2~ d )((d+l). This corresponds to a free fermion 
in any dimension with 

q+i(1, -1) = 1 - — , (free fermion) (117) 

For any n one also has the inversion relation: 

Lx n (-z) + (-) n Lr n (-l/z) = C n (118) 

where C n is a constant. For n odd, C n = 0, whereas for n even C n = —2(1 — 2 1 ~ n )((n). 
Using Lr n (0) = in the inversion relation, one finds: 

c rf+1 (oo,-l) = 2M -ij (dodd) (119) 

The above relation is fermionic, and since z = 00, it is an interacting theory with gap 
g = —00 and interaction parameter h = 1. As we explain in the next subsection, in 2D this 
theory is realized in e.g. the sinh-Gordon model with c qs = 1. In higher dimensions with d 
odd, c qs = 2(1 — 2~ d ) which is twice that of a free fermion and always more than a boson. 
We will refer to this the extreme-fermion since it possesses the highest possible c qs 

for a fermionic theory. 
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Notice that the above results were obtained without doing any integrals explicitly, but 
rather just by using functional relations satisfied by the Lr n (z). More complex polylogarith- 
mic ladders are constructed in the same spirit. Since the more recently discovered ladders 
are for N > 1, a classification of 1-particle theories is essentially within reach. In the next 
section we carry out this classification. 

X. RATIONAL THEORIES IN 2D 
A. One particle RSPE's 

Let us illustrate how ladders are constructed in the simple case of N = 1 (1-particle) and 
D = 2. The construction relies on additional functional relations which are valid at second 
order: 

Lv 2 (z) +Lr 2 (l-*) = C(2) 

Lr 2 (^) +Lr 2 (-z/(l - z)) = (120) 

Using the first relation above, one sees that Lr 2 (l/2) = £(2)/2. This, and the free fermion 
and free boson cases, were the relations known to Euler. Landen found more evaluations 
as follows. If r is a root to the polynomial equation z 2 + z = 1, then the above functional 
equations become linear equations for Lr 2 (z) with argument z — r, — r, — 1/r and r 2 . Let us 
take r = (\/5 — l)/2, which is the inverse of the golden ratio. This is the complete list of 
known RSPE's for D = 2 and N — 1, and has been known for over 220 years: 



Lr 2 (l) = 


-Lr 2 (- 


-oo) = 


= C(2) 


Lr 2 (l/2) = 


-Lr 2 (- 


-1) = 


i«2) 


Lr 2 (r 2 ) = 


-Lr 2 (- 


-r) = 




Lr 2 (r) = 


-Lr 2 (- 


-l/r) 


- 1^ 



Though it is generally believed that the above list is exhaustive, we could not find a proof 
in the literature. We suggest that a proof might be concocted based on proving that the 
known integrable scattering theories in 2D with one particle are exhaustive and are in 
correspondence with the above list. 
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Assuming that the relations eq. (jl21j) are indeed exhaustive, this leads to the classification 
of the possible theories shown in Table 1. In the table are listed the thermal gaps g and 
the interaction parameter h which can be inferred from the gap equation (|93|). The last 3 
columns will be explained below. 
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B. Integrable reconstruction in 2D 

Let us specialize the results of section V to integrable relativistic theories in 2D 3_5j . Let 
A a (9), a = 1, .., N denote particle creation operators where 9 is the rapidity parameterizing 
the one-particle relativistic dispersion relation: 

E = m a cosh 9 , k = m a sinh# (122) 

In this rapidity space the formulas for the partition function and integral gap equation take 
the form: 

/oo jn 
—m a cosh 9s a log(l - s a e-^ 0) ) (123) 
-oo 2vr 

i poo rn/ 

e a (6) = (macoshe - ti a ) + -J2 / — G ab (9 - 0> 6 log(l " s b e-^) (124) 

P ft J-oo 27T 

We will consider only theories where the scattering is diagonal. The factorizable S-matrix 
S a b(9) then enters the exchange relation as follows: 

A a (9)M0') = S ab (0 - e')A b {B')A a {B) (125) 
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In the case of 2D an alternative derivation of the kernel is possible [3J, |35(. Consider 
the quantum field theory on a space of 1 dimensional length L(= V). Requiring that the 
multi-particle wave-function is periodic leads to a quantization condition on the momenta 
k: 

e iLkl \[S ij {e i -e j ) = l (126) 

Taking the logarithm and differentiating with respect to the rapidity leads to a relation 
between pi and p as in eq. (|34|) with the kernel: 

G ab (9)=-i^\ogS ab (127) 

This derivation has its origins in the thermodynamic Bethe ansatz (TBA). In many cases 
one can construct explicitly the Bethe ansatz wave-functions and derive equations such as 
eq. (|12fi|) . For this reason we will refer to the above 2D equations as TBA equations. 



re last section, the kernel G can be deter- 
321 that if one defines the connected form 



On the other hand, using the arguments of t 
mined from the form factors. It was shown in 
factor in the proper way, then the resulting kernel is the same as in eq. ()127|) . 

We now attempt to identify the entries of the above table with specific integrable quantum 
field theories. All of the fermionic theories are well-known and their conventional names are 
in the last column. However each fermionic theory has a bosonic counterpart, which can 
be traced to the second functional relation eq. ()120|) . As explained in [35j], the statistics 
of the lagrangian fields does not necessarily correspond to the statistics in the partition 
function. First of all, there is no spin-statistics theorem in 2D. The S-matrix satisfies the 
unitarity condition S(9)S(-9) = 1. At 9 = 0, 5(0) 2 = 1 has two possibilities: 5(0) = ±1. 
The quantization condition eq. 1)126)1 imposes its own restrictions on the wave-functions. If 
5(0) = —1, two bosons are not allowed to have the same rapidity, and one should impose 
fermionic exclusion statistics. If 5(0) = 1, the situation is reversed. We can summarize this 
as follows. Let s/ = ±1 denote the statistics of the fields, where Sf — 1 is a boson. Let s 
denote the statistics of the partition function, as in eq. (j!23j) . Then s = s/5(0). Henceforth, 
in the 2D case, boson or fermion statistics will refer to the partition function statistics s. 

The reason that half the theories in Table 1 have not been identified is that nobody has 
ever made sense of an interacting bosonic TBA. It has even been suggested that interactions 
always turn bosons into fermions in 2D0. This issue was addressed in (4^, where a 
hypothetical model with S-matrix given by the well-known sinh-Gordon one, but with the 
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wrong statistics, was investigated. It was found that c qs was ill-defined in the UV, which 
shows one cannot simply change the statistics of known theories at will. On the other hand, 
as Table 1 shows, it is possible to have well-defined c qs for a bosonic TBA. 

For the remainer of this section we describe how to reconstruct S-matrices from the data 
in Table 1. We first review a few basic ingredients of diagonal factorizable S-matrices (see 



for instance 136 
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38j]). Unitarity and crossing symmetry require 

S ab (9)S ba (-9) = 1, S al (9) = S ab {m - 9) (128) 

where a is the charge conjugate to a. S-matrices satisfying the above equations are products 
of the fundamental building blocks: 

, sinh t;(9 + ina) 
fa{0) = . r . (129) 

sinh ^{9 — ma) 

If a particle is its own anti-particle, then crossing symmetry implies it is a product of the 
functions 

tanh ^(9 + ma) 

F a (9) = f a {9)U™ ~0) = : r fej ■ (130) 

tanh ^[9 — tira) 

The poles in the S-matrix are also severely constrained. If a simple pole exists in S ab at 
9 = iu, on the physical strip < iu < and with positive imaginary residue, then this 
corresponds to a stable bound state, with mass 

m 2 = m 2 a + ml + 2m a m b cos u (131) 

If this is a new particle, one must proceed to close the bootstrap. A negative residue pole is 
either the same particle in the crossed channel, or a sign of non-unitarity. We will need that 
f Q has a simple pole at 9 = ma with residue 2zsin7ra. From this it follows that F a has a 
simple pole at iira with residue 2itanna and a simple at i(l — a)TT with residue — 2itamra. 
The h values that appear in the gap equation for the Id integrable theories are: 

d9 d log S ab (9) 



i ab 



(132) 



-oo 27ri 39 

Using this one easily can compute that for each factor of f a , h — a — sign(a). For each F a , 
h = — sign(a). 

All the theories with h ^ are interacting. Let us first go over the identification of the 
interacting fermionic theories, all of which are known. 



32 



Fermionic theory with c qs = 1. Here since h = 1 the S-matrix must be _F_ a with a 
positive. For a positive, the pole on the physical strip has negative residue and thus does 
not correspond to a new particle, so the bootstrap is already closed. This is a one-parameter 
family of theories corresponding to the sinh-Gordon model, with lagrangian: 

S shG = J d 2 x Q(<90) 2 + A cosh^v 7 ^)^ (133) 

The S-matrix is known to be F_ a where a = b 2 /(2 + b 2 ). is a bosonic field, so Sf = 1, but 
since S(0) = —1, the partition function is fermionic. 

Fermionic theory with c qs = 2/5. Since here h = —1, the S-matrix must be F a with 
a positive. However for arbitrary a > 0, the physical pole will correspond to a new particle 
and this is no longer a one-particle theory. The only solution is to take a = 2/3. Then, 
since cos27r/3 = —1/2, the pole corresponds to a particle of the same mass as the original 
particle and can thus be identified with it. Thus the bootstrap closes with one particle. 
This model is known as the Lee- Yang model 48]. Generally, c qs = c — 24A m ; n where c is 
the Virasoro central charge and A m j n is the lowest conformal dimension of the fields. Here, 
c = —22/5 and A min = —1/5, giving c qs = 2/5. The residue of the pole is negative, which is 
an indication of the non-unitarity of the theory. 

Fermionic theory with c qs = 3/5. Here since h — 1/2 the S-matrix must be 

f-i/2 = — itanh \{6 — m/2). Though f a for general alpha cannot correspond to a 1-particle 
theory since it is not crossing symmetric, a = —1/2 is special since it is at least anti-crossing 



symmetric: f-iniin — 6) = —f_i/ 2 {9). This model is known 49| to be the $13 perturbation 



of the ^3,5 minimal model, with c = —3/5 and A min = —1/20, giving c qs = 3/5. 

We now turn to the bosonic theories. Note that they all occur at the same central charge 
as the fermionic ones. This does not imply they are the same theory, but only that they are 
the same at the UV fixed point. A theory is defined b y it s fixed point in the UV and by the 
operator which perturbs it away from the fixed point 50]. For the last 4 entries of Table 1, 
we only indicate the UV fixed point theory. 

Bosonic theory with c qs = 1/2. This is the most interesting case since the UV central 
charge is that of the Ising model, of which a great deal is already known. Let us attempt to 
reconstruct a consistent S-matrix for this theory deferring its interpretation. Since h = — 1, 
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the S-matrix must be F a with a positive. For the same reasons given above for Lee- Yang, 
a must be 2/3 otherwise the bootstrap doesn't close. We propose the following S-matrix: 

S(9)=pF 2/3 (6), p = ±l (134) 

We now carefully investigate the consistency of this theory. For p — — 1, the residue of 
the pole at 6 — 2ni/3 is positive, which is consistent with the theory being unitary. However 
the S-matrix now satisfies the bootstrap equation with an extra minus sign: 

S(6) = pS{9 + m/2>)S{9 - itt/3) (135) 

This minus sign can be removed by S — > — S, however this spoils the unitarity of the theory. 
For these reasons, the above S-matrix is not entirely consistent and its physical interpretation 
will remain an open question. The same is true for the next two bosonic theories. 
Bosonic theory with c qs = 2/5. Here h = — 2 which implies the S-matrix must be 
F a F a >. As in previous cases, the bootstrap will not close unless a = a' = 2/3: 

S(0) = F 2/3 (9) 2 (136) 

The resulting S-matrix has only double poles for the physical particle. Since 5(0) = 1, 
this model has a single bosonic field. Further analysis is needed to determine whether it is 
equivalent to the Lee- Yang model, or a different off-critical perturbation. 

Bosonic theory with c qs = 3/5. Since h = —1/2, the S-matrix must be: 

S(9) = -f 1/2 (6) (137) 

It has a pole at in/2, with negative residue, indicating non- unitarity. Here, Sf — 1. As for 
the fermionic version, this is an integrable perturbation of the non-unitary minimal 

model. Since the S-matrices are very similar in this case for the bosonic and fermionic 
theories (/1/2 = f-i/2(—@)), the two models are probably equivalent; again more analysis 
is needed. 

C. A supersymmetric 2 particle theory in 2D 

As was stated before, all of the newly discovered ladders are multi-particle. Not all of them 
can be brought into the RSPE form. It is beyond our scope to do a comprehensive study of 

34 



the vast number of results in so we concentrate on a simple example that doesn't appear 
to correspond to anything already known. Using the duplication and inversion relations on 
one of the results in 43], one can show 

c 2 (2/,l) + c 2 (j/,-l) = |, y = V2-l (138) 

This is an interesting relation since, unlike the known TBA's that are purely fermionic, this 
is a mixed theory, with one boson and one fermion, suggesting a supersymmetric theory. 

We can reconstruct the theory just from the data in the equation (|138j) . It's a two 
particle theory with z\ = z 2 = \/2 — 1, but with si = 1, s 2 = — 1. From the gap equation we 
determine the /i's to be: hu = h 22 = h\ 2 = h 2 \ = — 1. The choice of S-matrix: 

S ab (6) = -F 2/3 (6), Va,6=l,2 (139) 

yields these h values, and the bootstrap is closed for previously described reasons. The S- 
matrix shows no signs of non-unit arity. However the UV central charge is c = 3/4 which is 
not in the minimal unitary series. Our interpretation is that the physically sensible theory 
is two copies of the above. This has c = 3/2 with 4 particles and the supersymmetry 
is preserved. It is known that c = 3/2 is the highest central charge of minimal unitary 
super-conformal series jsij], supporting this idea. 

XI. RSPE'S IN 3D AND AD. 
A. One particle theories in 3D 

As explained in the last section, the second identity in eq. ()12U|) implies that in 2D a 
bosonic RSPE can be rewritten as a fermionic one with different gap parameter z. Because of 
bosonization, this is perhaps not surprising in 2D. In higher dimensions there is no analog of 
this identity so it is not possible in general to give both a bosonic and fermionic interpretation 
to a given RSPE. 



The only known 1-particle RSPE's (due to Landen 42j) beyond the free boson and free 



fermion in 3D have again been known for over 220 years and are the following: 

C,(il)4 c 3 (rM) = l (140) 
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where again r = (y/E — l)/2. This leads to the following table of possible 1-particle rational 
theories in 3D which is again very likely to be exhaustive, but we cannot prove it. 
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The interacting bosonic theory with c qs = 7/8 was suggested to be the 3D Ising model in 
section VIII. The c qs = 4/5 theory appears to be related to the O(N) sigma model, which 
at large N has c qs = f AfQ. The method for computing c qs in the latter work is just the 
effective action at large N, and is thus different than our methods. 



B. Multi-particle RSPE's in 3D 

For the remaining RSPE's presented in this paper we only reconstruct the particle spec- 
trum, gap equation, and interaction parameters h a b, since as previous stated we do not have 
the tools necessary to further reconstruct the theory. There are a number of generic features 
of these multi-particle RSPE's. They are generally of mixed statistics, i.e. involve both 
bosons and fermions. One can also check that in these examples, the interactions always 
reduce c qs from it's free field value, i.e. c qs < N b + (1 — l/2 d )Nf, where N b j are the num- 
bers of bosons and fermions. Below, we present this feature with the parameter r = 

^free 

However this is not generally the interactions in fermionic theories can increase c ns , 

the extreme-fermion being the limiting case. (See eq. ()162j) .) RSPE's appear to be more 
and more scarce as one moves up in dimension. This is not surprising since conformal field 
theories are also expected to be more scarce since there is no known underlying Virasoro 
algebra with a rich representation theory. 

Let us reconstruct the following c qs = y RSPE: 

13 

2c 3 (l/3,l) + c 3 (l/3,-l) = - (141) 
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It contains 2 bosonic particles with z\ = 1/3 and one fermionic particle with z 2 = 1/3. Here 
r = 26/33. The gap equations then read 

z 1 = (l-z l )- 2h ^(l + z 1 ) h ^, z 2 = {l + z 2 ) h ^(l-z 1 )- 2h ^ (142) 

The solution is h a b = —1, Va, b. 

Consider next the c qs = ^ RSPE with 3 bosons and 2 fermions with r = 11/18. 

3c 3 (x,l) + 2c 3 (x,-l) = ^, x = 3-2^ (143) 

Here the gap equation reads: 

z 1 = (l-z 1 )- 3hll (l + z 2 ) 2hl \ z 2 = (l + z 2 ) 2h ^{l-z 1 )-^ (144) 

The novel feature here is that there is a one-parameter family of solutions to h: 

7 7 u u 3/inlog(l - x) +logx 

ft 2 i = h u , h 12 = h 22 = — — — — r (145) 

21og(l + x) 

This suggests some exactly marginal directions since the fo's are required to be RG invariants. 
C. Multi-particle RSPE's in AD 

In 4D there are no known RSPE's beyond the free fermion, free boson and extreme- 
fermion. We present only two examples of multiparticle RSPE's. 
Here is an RSPE with c qs = 179/2 and r = 716/721: 

179 

42 c 4 (1/2, 1) + c 4 (1/8, -1) + 54 c 4 (2, -1) = — (146) 

It consists of 42 bosons and 55 fermions. Here there is a 3 parameter family of interaction 

parameters that give the gap parameters z: 

1 14 
h 12 = 14/in, h 21 = -y(h 22 - 1), h 13 = _ 27^11 

hi = -^(1 - 81h 33 ), h 23 = ~^h 22 , h 32 = -27/^33 (147) 

We leave the following c qs = 4393/8, r = 4393/4513 example as an exercise: 

4393 

264c 4 (x, 1) + 15c 4 (x 2 , -1) + 4c 4 (x 3 , -1) + 324c 4 (l/ X , -1) = — (148) 
where \ = 2 — \/3. 
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XII. FRACTIONAL STATISTICS AND IMAGINARY CHEMICAL POTENTIAL 



A. General arguments 

In the formalism we have developed so far, the bosonic or fermionic statistics (s = ±1) 
has played a crucial role. It is well-known that in 2D and 3D, fractional statistics particles, 
i.e. anyons, are also possible[53]. In this section we argue that this possibility can be 
incorporated with an imaginary chemical potential. The connection between anyons and 
an imaginary chemical potential has previously been suggested in the context of the 3D 
Gross- Neveu model (3). 

Let us model our discussion based on the 3D case, though we will also present results in 
2D. The standard way to obtain anyons is to add to the action a topological, or Hopf, term: 

S^S + i$S top (149) 

The simplest case, which we will refer to as the free anyon, corresponds to S being free 
charged bosons or fermions with U(l) current J M and charges Q — ±1. Here S tD p is the 
Chern-Simons term where the U(l) gauge field is manufactured from the U(l) matter 
current = e flua d u A a . The statistical parameter $ is normalized so that the spin of the 
particle is 19/271. 

At finite temperature, space-time has the topology of R d <S> S 1 , where S 1 is the one 
dimensional circle. Let C denote the configuration space of the field theory; for example for 
the 0(3) sigma model, C = S 2 . At finite temperature the relevant homotopy group is the 
one classifying maps from R d <g> S* 1 — > C which we will denote as ttt- 

Let us assume ttt = Z and the Hopf term S top is ±1 for the particles of U(l) charge 
Q = ±1. This implies that the particles contribute ±ii9 to logZ. Comparing with eq. (|35j) . 
this corresponds to a chemical potential /i = ±ii?T where T is the temperature. More 
generally, let \x a denote the imaginary chemical potential which modifies the statistics of the 
a-th particle with statistical parameter 9 a : 

Va = l<&aT (150) 

The important feature of this statistical chemical potential is that \i a jT remains finite at 
any temperature, which is a consequence of the fact that S top is finite. This leads to a 
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finite contribution to the thermal gaps g a = —i$ a - In terms of the gap variables z a the gap 
equation reads: 

z a = - s b z b )- h ^» (151) 

b 

where as before s a — ±1 corresponds to whether the original undressed particles were 
fermions or bosons. 

If the particles are free, h ab = and z a = e i1?a , i.e. the thermal gaps are purely statistical. 
Let us denote by cd{$, s) the value of c qs in D-dimensions where s — 1 (s = — 1) corresponds 
to ^-statistically dressed charged free bosons (fermions) with $± = ±-#. From the results of 
section VII, we propose the following formula for these "free anyons" : 

cd(^, s) = ^ (Lr D (se il? ) + hr D {se^)) (152) 

(For z on the unit circle, \z\ = 1 in the formula ()103|) so that Lr(^) = Li(z).) 



B. Free anyons in 2D 

A non-trivial check of eq. (|152j) is available in 2D due to the remarkable dilogarithm 
identities for z on the unit circle: 

c 2 (tf,l) =2 (1 - — (1 - ^- ) ) , 0<^<2tt (153) 



7T \ 27T 



c 2 (i?,-l) = 1-12 ( — ] -7r<$<7r (154) 



The arguments leading to eq. ()152j) assumed the underlying field theory was a free 
theory of charged fractional statistics particles. In 2D the spin a of a field can be defined 
by monodromy property: 

^(e 2 ^z, e~ 2m z) = e 2 ™i){z, z) (155) 

where z,z are the euclidean light-cone space-time variables, z = x + w,z = x — iy. The 
simplest free theories of anyons in 2D are the so-called (3 — 7 systems [54] which were first 
studied in connection with the world-sheet ghosts in string theory. The theories are defined 
by the action 

S = 27 / (P 9 * 7 + P 9z ^ (156) 
where the spin of 7) is (cr, 1 — a) and the opposite sign for (/3, 7). The fields can be 
quantized as bosons or fermions, the latter usually being referred to as a b — c system. 
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Consider first the bosonic (3 — 7 system. The central charge is known to be 3| : 

c ;3 _ 7 = 2(6(T 2 -6(7 + 1) (157) 

Due to the identity eq. ()153|) . this corresponds precisely to the statistically dressed boson, 
i.e. C/3_ 7 = c 2 ($, 1) with $/2n = a. 

Consider next the case where (5 — 7 are quantized as fermions, usually denoted b — c. The 
central charge is known to be: 

c b _ c = -2(6a 2 -6(7 + 1) (158) 

Since the undressed particles are already fermions, by the same reasoning as before, this 
should correspond to a statistically dressed fermion with $/27r = a — 1/2; it is easily con- 
firmed that 02(1?, — 1) = Cb-c using eq. ()154j) . 

It is somewhat remarkable that we have been able to obtain the above central charges 
simply from a statistical thermal gap. We point out that equality of cg_ 7 , q,_ c with c 2 ($, s) 
is only valid for the range of $ indicated in eqs. ()153I154|) since 02(1?, s) is a periodic function 
of $ whereas cg_ 7 , q,_ c are quadratic. 



C. Rational free anyons in 3D 

As the equations ()1531 1154J) show, when $/27r is rational then 02(1?, s) is rational. In 3D, 
there is no analog of these equations and c^i?, s) is generally irrational. However for some 
special values of C3 turns out to be rational. 

Anyons of spin 1/3 or 1/6 are examples of rational theories: 

< »(3- 1 )=-«»(t-- 1 J = 3' ^(s'-^-^iT'V^S (159) 
An interesting feature of this case is that these gap parameters can be obtained from an 
interacting theory with no statistical chemical potential. Consider for example the bosonic 
theory with C3 = 2/3. There are two particles with gap parameters z% = e ln ^ 3 ,Z2 = e~ l7T ^ 3 . 
These are the solution to the bosonic gap equation 

z x = (1 - Zl )- h "(l - z 2 )- h ^\ z 2 = (l- z x y h ^{\ - z 2 y h ™ (160) 

with real interaction parameters h\\ = h 22 = 1/2, hi 2 = h 2 i = —1/2. This observation 
suggests that anyons can arise in an ordinary interacting theory without inclusion of a Hopf 
term. 
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Another example we have found is for particles of spin 1/4: 

* (*_!)=_*('!) =3. (161 ) 

It remains unknown to us whether there are other values of d that correspond to rational 
theories. 



XIII. PHASE STRUCTURE 

In this section we reintroduce the chemical potential and address a basic issue: for given 
interaction parameter h and chemical potential fi, is there a fixed point with a well-defined 
c qs ? Generally speaking, /x is related to the density of particles, where \x > (/i < 0) 
corresponds to increased (decreased) density. This issue has largely been ignored even in 
the 2D integrable TBA approach. 

Given the interaction parameters h and the solutions to the gap equation z, c qs as given 
in eq. (jl(J5|) is not necessarily real, or even well-defined, and this could signify something 
like a phase transition. In this section we explore this issue for 1-particle theories, where a 
real c qs requires a real z. Of course as the last section shows, real c qs can arise for complex 
z in multi-particle theories. 

The functions Lr^ +1 (^) are real for — oo < z < 1. This implies the following physical 
ranges of z: 

bosons : — oc < ; < 1 . 2 ( 1 — — 1 < c qs (z) < I 



2 d 



fermions : - 1 < z < oo, - 1 < c qs (z) < 2 M. - (162) 



In both cases c qs < corresponds to z < 0. 
Let us write the gap equation as 



where /i is the chemical potential. We are only considering fixed points here, so fi must be 
proportional to T for it to have any affect on c qs . Given h, //, one can graphically obtain the 
gap z by first plotting the function f2 (z) as a function of z for the ranges in eq. (j!62)l . A 
given /i, by eq. (|163jh is a point on the y-axis. Thus, the intersection of a horizontal line 
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crossing the y-axis at e^/ T with the curves determines the gap z. There are several cases 
depending on h and whether the particles are bosons or fermions. 

Bosons with h < 0. This case is shown in figure 2. One sees that for any real \x there is 
a single solution with < z < 1. Solutions with a negative z, on the negative y-axis, require 
the chemical potential to have an imaginary part Im(fi) = nT. As explained in the last 
section, this can correspond to a statistical chemical potential with i? = tt. In this regime 
there are two solutions to z if h < — 1, one solution if — 1 < h < 0, and no solution if the 
real part of \i is too large. Since in all cases a negative z requires an imaginary part to the 
chemical potential, throughout this section by "complex chemical potential" we mean one 
where specifically Im{ji) = nT. 











1 

l+h 










1 



FIG. 2: Bosons with h<0. 

Bosons with h > 0. This case is shown in figure 3. Here, for real chemical potential, 
there is no solution unless fi is negative and 

H/T< - ((/i + 1) log(/i + 1) -hlogh) (164) 

in which case there are two solutions with < z < 1. Thus there is no solution for 
zero chemical potential. The horizontal line displaying the two solutions is rendered in the 
figure. This implies that if the gas is too dense, there is no fixed point. The latter means 
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for instance that the black body formula breaks down. This could signify a kind of Bose- 
Einstein condensation. This explains for example why one cannot obtain a well-defined UV 
limit for the bosonic version of the sinh-Gordon model, as explored in j^, unless one lowers 
the density. As the density is decreased, there are two possible c qs . It is not entirely clear 
what the significance is of these two solutions. If they are both possible in the same model, 
then this suggests there are two phases and which one is reached depends on the details 
of the RG flow. On the other hand there is always a solution with z < with a complex 
chemical potential. 




z 



FIG. 3: Bosons with h > 0. 

Fermions with h < 0. See figure 4. For any real \x there is a single solution with 
< z < oo. With complex chemical potential there are two solutions if Re(fi) is not too 
large. 

Fermions with h > 0. See figure 5. For < h < 1, and real negative chemical potential 
there is one solution with < z < oo. The sinh-Gordon model is the special case h = 1 
where there is one solution with z = oo. For h > 1 there is no solution unless 

H/T < -(hlogh- (h- l)log(7i-l)) (165) 

in which case there are two solutions. Finally with complex chemical potential there is 
always a unique solution with — 1 < z < 0. 
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FIG. 4: Fermions with h < 0. 
XIV. CLOSING REMARKS 



The quasi-particle form of the finite temperature partition function and the integral gap 
equation is a new formulation of finite temperature field theory that appears to be well 
suited for capturing certain phenomenon in a clear fashion. Though the method is an 
approximation in general, and the nature of this approximation needs further clarification, 
we believe that it captures the essential mechanism that leads to non-trivial fixed points: 
interactions lead to thermal gaps that are proportional to the temperature, A = gT, and c qs 
is given in terms of some very special combinations of polylogarithmic functions eq. of 
the complex variable z = e~ 9 . The spectacular properties of these functions Lr n (z) lead us 
to believe that this is a general property, i.e. that any non-trivial fixed point has c qs given 
by the expression f)109j) . The classification of rational conformal theories in any dimension 
is then essentially mapped onto the problem of the construction of polylogarithmic ladders 
in mathematics. It would be very interesting to use this idea to explain some of the rational 
values of c qs obtained for super conformal gauge theories using string techniques 

The thermal gaps are related to the fundamental interaction parameters of the theory by 



44 



e /VT 




-1 






Z 






h-l 













FIG. 5: Fermions with h > 0. 

the algebraic gap equation (|1()8|1 . There may be other methods that lead to the algebraic 
gap equation besides the integral gap equation proposed in section V, such as the finite size 
effective potential I7|, and this is important to explore. 

One of the main questions left unanswered in this paper is how the structure of the stress- 
tensor form factors can lead to finite interaction parameters h^b that are RG invariants. In 
2D this is well-understood using integrability and the relation of the form factors to the exact 
S-matrix. Understanding this is clearly necessary for further progress on specific models. 

We have touched only briefly on the possible applications. Further work is clearly needed 
on some of the suggestions made in this paper. For example, the proposals for the 3D Ising 
model and the free energy of anyons in 3D should be compared with numerical simulations 
of the free energy. 

To phrase properties of the equation of state of matter in cosmology in terms of the 
c-function C qs (T) appears to be fruitful. As we argued in section IV, this leads to the 
bounds —1 < w < —1/3 for dark energy, but perhaps more importantly we have related 
the temperature dependence of w to the RG flow of C qs , and thus "c-theorems" become 
relevant. For instance, from this perspective the fate of the accelerated expansion depends 
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on whether the dark energy has an infra-red fixed point. 
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